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Abstract 

High derivative terms do not play a major role in field theories because of the 
associated complexity and inherent difficulty in connecting these terms to phys- 
ically measurable quantities. A role for higher derivative terms is analyzed for 
the case of field theories used to describe phase separated systems. In these the- 
ories, higher derivative terms are directly connected to an interfacial free energy 
which contains the mean and the Gaussian curvature and are shown to determine 
explicitly the shape of the interface. 
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1 Introduction 



Cooperative behavior of lipid molecules in biological membranes is believed to play a cen- 
tral role in processes happening at cellular level. More generally, cooperative processes 
appear in systems exhibiting hydrophobic-hydrophilic interaction even beyond those ap- 
pearing in biology H]. Cooperative processes are captured and studied with difficulty 
by numerical simulations at molecular level, mainly because those methods have a low 
probability of capturing events involving multiple correlated molecular steps. Moreover, 
simulations at molecular level are usually plagued by an extensive range of parameters 
and although successful in predicting short-time molecular dynamics, they make less of 
an impact in a realistic statistical system that covers large spatial and temporal scales. 
Other computational approaches are bound to improve the relatively straightforward and 
parameter-extensive all-atom simulations. These methods made their way towards molec- 
ular simulations [2] and they appear under the terminology of coarse-grained numerical 
methods. Coarse-grained methods are based on a considerably reduced parameter space 
and provide a significantly increased speed. Therefore, coarse-grained methods cover the 
long-time dynamics familiar to many cell-biology cooperative processes (nanoseconds to 
seconds) with the additional advantage of allowing an intrinsic inclusion of statistical aver- 
aging. In spite of the advantages, coarse-grained methods present challenges of their own. 
For example, in the framework of these methods, making contact with realistic systems 
and identifying those few most relevant parameters is usually a difficult task. 

Particular chemical systems where cooperative behavior is responsible for a geometri- 
cal rich phase structure are polymeric blends. For this case, the cooperative behavior is due 
to the long-range interaction stemming from "stretching" abilities of polymer chains. A 
coarse-grain computational and theoretical model where molecular details and coordinates 
are replaced by a field theory of concentrations is already applied and used for a while 15], 
lU . Self-consistent field theory methods which are based on a mean field approach provide 
a successful representation of the experimental data. These methods were recently applied 
with success to biological membranes and other complex cellular structures, such as vesi- 
cle, see for a review fS). One of the major drawbacks of the "self-consistent" methodology 
is the lack of transparency. In particular, a direct relation to physical measurable param- 
eters, such as interfacial tension and geometrical invariants, such as mean and Gaussian 
curvature of interfaces is difficult to write explicitly. Moreover, the method faces diffi- 
culties when extended to dynamical models 0, difficulties similar in magnitude to those 
appearing in molecular models. 

A more explicit field theory model can be obtained for a system of homopolymers, us- 
ing a formal perturbation expansion of the corresponding self-consistent field theory. The 
result of the expansion is a regular free energy expression that contains higher derivative 
free energy terms Q . This type of terms come regularly in effective field theories from 
integration of higher energy modes (or alternatively of quickly relaxing modes). Although 
these terms appear in other field theories used to describe strongly correlated systems |8 1, 
they are not regularly treated mainly because the complexity and lack of transparency as- 
sociated with them. It is the goal of this work to clarify the connection between higher 
derivative free energy terms for a phase separated system and the coefficients that govern 
the dynamics of interface geometry. 
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A common phenomenological form for the interfacial action was proposed by Helfrich 



where 7 is the interfacial tension, ch is the free energy coefficient associated with mean 
curvature H, Hq is the intrinsic mean curvature and cq is the free energy coefficient as- 
sociated with Gaussian curvature K. Such an interfacial free energy expression is useful 
both from the experimental and theoretical point of view, but the above free energy can 
be used to describe a static interface and does not allow a continuous representation for 
common experimental processes such as joining and splitting of interfaces. In this study, 
starting from a generic free energy that contains additional higher derivative terms, the di- 
rect contact with the free energy ([T]i is made. In addition, it is shown that the interfacial free 
energy naturally features contributions from the interfacial width, contributions that can be 
directly correlated with the disappearance or formation of interfaces. Another advantage 
is that from the point of view of the number of degrees of freedom, a complete match- 
ing between bulk and interfacial degrees of freedom is achieved. Matching of degrees of 
freedom is reached through the usual scheme of dimensional transmutation, where a coor- 
dinate, the one that is normal on the interface, is transformed in an additional scalar field 
that lives on the interface, the interfacial width. The picture is not yet complete, but the in- 
terfacial width presents itself as the required intermediary field that provides a continuous 
description to the nontrivial topological change associated with joining and splitting. The 
interfacial free energy thus obtained reveals also nontrivial interactions between the inter- 
facial width, and the geometrical properties of the interface, such as mean and Gaussian 
curvature. Such interaction terms come from higher derivative terms in free energy and 
open the possibility of modulating the interface itself, in particular the interfacial width, 
by varying the curvature of the interface. 

The plan of the paper is as follow. In section 2, the method of derivation for the inter- 
facial free energy from a generic free energy action is presented. In section 3, the method 
is applied to the case of Landau-Ginzburg 0^ potential and all interfacial coefficients are 
calculated. The general case is presented in section 4 and section 5 contains the calcu- 
lation of interfacial parameters for the realistic case of binary incompressible blend of 
homopolymers with different molecular weights. 

2 From higher derivatives to interfacial properties 

Higher derivative terms appear naturally in effective free energy actions for statistical sys- 
tems and their appearance stems from an integration of degrees of freedom with relaxation 
time scales shorter than the scale time of interest (namely, those degree of freedom which 
are inaccessible to the particular experimental setup). Alternatively, these terms can also 
appear in quantum effective actions, but here, their appearance is a sign of the underlying 
intrinsic quantum dynamics of fields. The problem becomes more involved when the free 
energy admits soliton-type solutions, solutions that in a traditional physics setup excite 
both low and high energy modes and it is not always clear how to separate (and integrate) 
the energy modes. 



El: 




(1) 
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In this work, it is assumed that following an integration of inaccessible modes, an ef- 
fective free energy with higher derivative terms is already available. This effective free en- 
ergy still possesses an effective potential that admits domain - wall solutions (soliton-like 
solutions). Within this setup, a free energy admitting domain-wall solutions and having 
also higher derivative terms, the impact of higher derivatives terms on the physical char- 
acteristics of the domain-wall solution is analyzed. Although such a free energy may look 
esoteric, it is shown in section 5 that such expressions can be obtained in realistic systems 
such as an incompressible mixture of homopolymers. 

The generic free energy expression considered in this work is: 

= F^[4>]+F2[4>] (2) 
Fi[4>] - j d\{v{c^)+M{c^){Vd^f) (3) 



+RmW^<l>)^ + 5(0)9„(/.c)'J'(/.) (4) 

where Fi in equation (O, is the low energy part of free energy and F2, in equation (|4]i is 
the correction part of the free energy coming from integration of degrees of freedom and 
contains higher derivatives. F2 is considered to formally be a perturbation effect added to 
Fi. V{(j)) is a positive potential that has two minima ^'_,+ corresponding to the two phase 
components of the system. The value of the potential in those phases, namely at minima 
points ^'-,+, is considered to be zero. Under this assumption, the free energy density in 
the bulk of each phase, namely in the regions where the field is constant and equal with 
either or '5+, is 0. Thus, no contributions to the free energy come from bulk regions 
and in fact, all free energy contributions come from the spatial region where the values of 
the field varies between the two coexistence values ^'-,+. These suggest that a complete 
reduction of free energy to such a region can be made giving rise to a free energy defined 
on a one dimensional lower space, namely defined on the interface. In order, to employ 
such a reduction, an Ansatz using the domain-wall solution of the equation of motion for 
the first two terms of the free energy, Fi is used. This equation of motion is: 

- 2Af(0)-^0 - M'{4>){^(pf + V'{4>) = (5) 
ax'^ ax 

where the ' denotes the first order derivative in (f). Following the usual procedure for 
domain-walls, this equation can be reduced to a first order differential equation: 

Tx'^^y-m) 

The solution to this equation that describes a domain wall with an interface at x = can 
in principle be determined and denoted as '^{x). This solution goes from ^ at x = ~oo 
to 5*+ at X = +00. The maximum of the first derivative of '^{x) is taken to be at a; = 
and this point is considered to be the interface location. The following field redefinition is 
proposed: 

4>{r) = vE'(/(r)) (7) 
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where /(r) is a field with dimension of length. The field /(r) represents physically the 
geometric profile of the interface, and it is used to describe the location of the interface 
in three dimensions through the equation /(r) = 0. The free energy is rewritten in terms 
of the geometrical field /(r) used in O using the relations below between derivatives of 
field (f) and those of field /: 

d^<P = ^'{f{r))d^f{r) (8) 
dc.p<t> = ^'{f{r))da,pf{r) + ^"{f{r))d„f{r)dpf{r) (9) 

These equations are substituted in equation (|2]i. A simplification in calculations is brought 
by expressions that relate the derivatives of '!''(/) and ^'"(/) in terms of ^'(/) obtained 
using the equation (|6]l and equations derived by its differentiation: 

- (^y (10) 



^ 1 f V'{^)M{^)-V{^)M'{^) \ 

2\ M(*)2 J ^ ' 

where depends on / and represents the domain-wall profile. In terms of the field /(r) 
and its derivatives the free energy can be finally written as: 

F[f] = F,[f] + F,[f] (12) 
=jd\M(^)^'^(l + {Vff) 

F2[f] =J rfV(^(7V(*)*'4 + (P(*) + Q(*))*'2^" + (i?(^) + 5(*))*"2)(V/)4 + 



It is worth mentioning that in the equation above, the dependence on / comes both explic- 
itly, as seen above, and implicitly through the profile function ^(/). The free energy in 
terms of the field / are similar in terms of derivative term representation with the initial 
free energy. Further simplifications are obtained by choosing a particular Ansatz for /. 
The choice in / removes all geometrical ambiguities related with describing the surface 
through a functional equation of type /(r) =0 and allows the introduction of interfacial 
physical parameters in the free energy. A local system of coordinates is used to simplify 
the calculations of interfacial free energy terms. The derivation is done for a particular in- 
terfacial point and a particular system of coordinates. The expected geometrical invariance 
of interfacial free energy is used then to generalize the local free energy to a geometrically 
invariant form, valid at any interfacial point and in any coordinate system. The analysis is 
carried in three dimensions, but it can be observed that most of the features can be gen- 
eralized directly to any dimension. By translation and rotation transformations, the given 
interfacial point is considered to be (0, 0, 0) and the Ansatz for f{x,y, z) around this point 
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is: 

f{x,y,z) = X{y,z)(^x ~ a{y,z)^ (13) 

cr(0,0) = 
9j,,,(t(0,0) = 
A(0,0) > 

where x — (T{y, z) is the location of the interface in the new parameterization and \{y, z) is 
a scalar field on the interface that can take only positive values and is directly proportional 
with the inverse of the interfacial width. Such a choice for /(x, y, z) although apparently 
limiting, in fact, does not make any restriction on the interface. Therefore, the choice in 
(fT3T l represents the most general interface in a particular local coordinate system valid in 
the neighborhood of the point (0, 0, 0). The choice can be justified by trying to fix the 
large freedom in choosing the curve equation as a function of regular spatial coordinates. 
For example, the interfacial function / is in first approximation linearly proportional to 
the normal coordinate x at the interfacial point (0,0). The final form for / in equation 
(fTST i is obtained through translation and rotation of the local coordinate system such that 
X axis is normal on the surface x = a{y, z), whereas y, z axes are parallel to the tangent 
plane of interface at interfacial point (0,0). The free energy density is computed in the 
neighborhood and the limit {y, z) — > (0, 0) is used freely in the computation. Overall, the 
procedure outlined above amounts to reduction of the free energy density along the direc- 
tion normal to the interface using the known expression of the one dimensional profile. As 
it is shown in the appendix 1, at point (0, 0) only, the metric in parameterization ( fTsT l. is 
the trivial two dimensional metric. The second fundamental form is given in terms of sec- 
ond order derivatives of a, whereas the mean curvature is given by V^cr. In order to fully 
reduce the three dimensional free energy to a two dimensional interfacial free energy, it 
is necessary to absorb the normal direction x. This is accomplished by integration of free 
energy density at each point of the interface, on normal coordinate to the interface, x. The 
coordinates y, z describe the interface and the assumption is made that after integration a 
local interfacial density is obtained. A straightforward integration over the normal coordi- 
nate, X is done by a change of variables in a neighborhood of the point (0, O)from x to /: 
X = jj^j-^ + (^{y, z). At the interfacial point, (t(0, 0) = and x = jjjj-q^- This change 
of coordinates leads to the following expression for integration ((y, z) — > (0, 0)): 

I dxdydz... — > f dydz ( , ... (14) 

J J J Ky,z) 

Integration over J df is made first at the point (y, z) — (0, 0) and it is followed by inte- 
gration over the interface coordinates y and z. As it was mentioned above, the integrand 
expression obtained at (0, 0) in the chosen local coordinate system is generalized using the 
expected geometrical invariance of the final expression of free energy density. Thus, the 
trivial metric and measure at point (0, 0) becomes the full metric on the interface, g with 
gaf3 ( a and f3 are interfacial coordinates, e.g. y and z) and measure c?S = ^/detgdydz. 
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For completeness, the expressions obtained for each term in equation (|2|l are given 
below: 



(VA)2n2 



A2 



+ 2i?(*)*'*") f A' + f^^^) f/^ - 2\H 



iV(*)*'4 + ( + Q{-^) ] *'2*" + [R{^,) + 5'(^-) ) ) ( a2 + ) + 

7A)2w 
A2 A-' A 

/ V2A \2 

+5(vi/)*'2(/2^^£^f!!^ + 2(VA)2 - 2/a„^An"'5 + A^n^^n"^)) (i5) 

where 5* depending on / is the domain-wall solution in the absence of higher derivative 
corrections, 5*' is its first derivative and 'I'" is its second order derivative (see equation 
©, ( flOb . (fTTTi). A depending only on interfacial coordinates, e.g. y and z, is the inverse 
of interfacial depth at the interfacial point (y, z). 11"'^ is the second fundamental form of 
the surface as embedded in the three dimensional space and is related with mean curvature 
H = ^tr{Pi) and Gaussian curvature K = dct(n) through equations: 

H = (16) 

m'^~2K ^ no^pTP^^ (17) 

(see also appendix 1). In the above expression, all derivatives are understood in terms 
of geometrical derivatives in interfacial coordinates, namely da- The dependence of A 
and interfacial geometry is completely featured and only the corresponding values for 
coefficients are still to be determined by integration over /. The dependence on / in the 
integrand comes directly and indirectly through the functions '!'(/), ^''(/) and 
Solving for the exact contribution to interfacial terms from each term implies integrating 
over / and that can be done in the case when the exact expression of interfacial ^(/) is 
known, namely for simple forms of potential V{(t)), see section 3. Moreover, calculating 
the coefficients can also be made in other particular cases, when the potential is more 
complex and an exact analytic expression for is not possible, see section 5. An 

interesting feature of equation (fTSI l is that it pinpoints the contribution of each of the terms 
from dU to the interfacial free energy ([T] The summary of these dependencies is shown 
below for each geometrical term appearing in (flSl l 

M{4,), N{4,), P{4>), Q{4>), R{ct>), S{^)^X 

R{(j}) — > 



1 



where the first expression shows the contributions to terms containing only A, whereas 
the other expressions show contributions to terms proportional with H, and so on. 
The terms multiplied by M {</>) and N{(f)) do not bring curvature terms in the interfacial 
free energy. Terms proportional with K come only from S'((/)) term, with from S{(t)) 
and terms, whereas the intrinsic curvature Hq is modulated by Q{4>), R{4>) 

and S{(f). The expression (fTSl l shows that in addition to the geometry of the interface, 
the free energy density also depends on a scalar field, the inverse of interfacial width A. 
The geometric invariants that make their way in the free energy explicitly at this level of 
approximation are those enclosed in the first fundamental form (metric) and second fun- 
damental form, namely 11"'', or indirectly in the better known invariants, mean curvature 
H and Gaussian curvature K. Since the goal of this work is to make contact with the in- 
terfacial free energy described in equation ([T]), only the terms of order are considered 
(|2]l. It is worth mentioning that the same method can be applied for even higher derivative 
terms, e.g. of order V^, resulting in the appearance of other geometrical invariants in the 
interfacial free energy. 



3 Simplified example - Symmetrical double-well potential 

The expression developed in the previous section is applied to a simplified case. In this 
case a Landau-Ginzburg type free energy with a symmetric double-well potential, V{4)), is 
considered and higher derivative terms are added. Furthermore, for this section coefficients 
in the free energy equation (|2]l are constant functions. In this case, the terms that come 
multiplied by Q and S are in fact redundant, their contributions being absorbed, up to 
total derivatives, in P and R constants. Also for clarity, the one dimensional domain-wall 
solution considered as reference in equation (|7]l is taken as simple as possible; 

*(a;) = tanh(x) (18) 

This expression for the domain-wall solution is consistent with values of the coefficients 
in M and V{(t>) shown below. The free energy expression considered in this section is: 

= \{cf'~lf; A/(0) = i; ^ = (19) 



F[<P] = J _ 1)2 + 1 (V0)2 + iv(V# + P(V</.)2(V2^) + R{V 

Equations ( fTSl ) and the explicit form for domain-wall solution, equation ( fTSl l allow the 
explicit integration along the normal direction over /. The full derivation is shown in the 
appendix 2 and only final results are given below. The first two terms in the free energy 
expression ( fT9] l lead after field redefinition, geometric identifications and integration over 
/to: 

'2 1 + A2 7r2-6(VA)2 



Jdydzd^ 



.3 A 18 A3 

The term multiplied by coefficient N leads to: 

'32 3 8(67r2 - 49) (VA)2 2(37r4 - SStt^ + 60) (VA)"^ 



(20) 



N 



dydz[^-\ + — —+ — —j (21) 
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and the term multiplied by coefficient P shows the following contribution: 

32 2 2{30 - 47r2) (VA)2> 



P 



^yH--5'+ 45 )^ ^"'^ 



Finally, for the term whose proportionality constant is R, after integration over / and 
cancellation of total derivative terms, the following expression is obtained: 



15 45 A 225 A^ 

A3 — + — — A^ + y^^J ^^^^ 

After collecting all terms from equations ( |20l ). ( ISTT i. (l22l i and (|23])and generalizing the 
expression from the neighborhood of interfacial point (0, 0) to any point on the interface, 
the following expression is obtained: 



(VA)V7r^-6 ^8(67r2-49) 8tt 
^ ' ' ' N— - + R 



■(■ 



A V 18A2 315 45 

(VA)4/ 2(37r4- 357r2 + 60) 771^-360 

+ - — ^{N— + R 

A5 V 225 225 

6-7r^ (VA)^V^A ^^-6 (V^A)^ 

A3 A^ + 

The above explicit interfacial free energy expression shows the nontrivial character 
of the higher derivative terms for the case of a simplified potential. It shows a rational 
dependence on interfacial width, which can not be obtained through symmetry or other 
qualitative arguments. Furthermore, the free energy for this particular case does not de- 
pend on Gaussian curvature, but only on mean curvature. The independence of free energy 
on Gaussian curvature is consistent with the fact that S is constant function and with the 
topological character of two dimensional Gaussian curvature. Gaussian curvature appears 
only when a more involved function is chosen as S{4)) in equation (|2|l. Lastly, the expres- 
sion (l24b shows explicitly the coupling between the geometric parameters such as mean 
curvature and interfacial width suggesting that it is possible to affect interfacial width, e.g. 
extend or contract, by manipulating the curvature of the interface. More detailed conclu- 
sions will be given in the next section and in section 5, in the case of homopolymers blend 
where entropic free energy contributions for internal degrees of freedom for an extended 
object are responsible for the higher derivative terms. 



4 Interfacial free energy - General case 

To gain a better understanding for the interfacial free energy, a few assumptions are made. 
The expression (fTsT i is simplified to a level that allows an analysis of the effects of higher 
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derivative terms. The main assumption is that the interfacial width does not vary across 
the interface, meaning that VA = and V^A = and that results in the elimination of 
such terms from free energy expression ( fTSl l. This assumption restricts the type of systems 
where the method can be applied. The assumption is also satisfied for a general system, 
but in cases when the interface has an intrinsic symmetry. For example, planar, spherical, 
cylindrical or toroidal interfaces in three dimensions have such symmetries. Under these 
conditions, the free energy expression becomes: 

F = Jd^jJ d/(^M(1')^''2(^l + A2) + (iV(1')^''4+ (p(*) + Q(«'))^''2*" + 

+ 5'(*))l'"2) A"* - 2(^P(^')*'3 + 2R{^)^'^"^X^H + 

+4:R{^)'i''^X^H^ +'IS{^)-^'^X^(^H^ - (25) 

where equations ( fT6] l. ( fTTI ) are used to provide transparency. The assumption taken at 
the beginning of the section leads to a free energy density that does not depend on / 
explicitly, but through It is possible now to avoid solving the domain-wall equation 

^ explicitly, task which proves to be difficult most of the times. The following change of 
coordinates is made: 

/ ^ * = *(/) (26) 

where is the domain wall solution of the field theory and thus / is given in terms of 
the inverse function of ^(/). This leads to the following change for integration variables: 

d.dyd^... J dyd^ J j§;^y ^ J § - (27) 

where is expressed in terms of in equations ( flOb . dTTT i. Using equation dZSl ) and 

the limiting procedure dZTl i results in: 



F[X; <j] = j d^ (^Ai (a-1 + a) + AaA-'' 



+A^X^(^-H)+AiXH^ + A^x(H^-^^^ (28) 

where a{y, z) is a parameterization of the interface and A(y, z) is the inverse of the inter- 
facial width in the normal direction, and dS is the geometric measure of integration. The 
expression for coefficients after the change of variable to '5, .5, are: 

/•*+ 

Ai ^ j d«'M(*)*' 

A3 = 2 y d'^{p{^)-^''^ + 2R{^)-^"^ (29) 
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/■*+ 



+ 



where Af, N, P, Q, R and S are the functions of (/> in the original free energy ^ and the 
and 4*" denote the first order and second order derivatives of the in terms of ^I^. 
All these coefficients are defined with the normal to the surface oriented from the phase 
having to the one with 5'+. The identification of parameters using equation ([T]i and 
(|28] | gives; 

7 = Ai(a-i + a) + AaA^ (30) 
CH = (A4 + A5)a (31) 
CG = -^A (32) 

^ ^(i^h^^ ^^^^ 

These expressions connect the free energy coefficients generically present in the equa- 
tion (|2]i and the interfacial expression proposed phenomenologically in When the 
additional derivative terms are absent, the free energy reduces to the first term of Ai leav- 
ing the second term of Ai and A2. .5 zero. Only a contribution from interfacial area and 
only the interfacial tension parameter is nonzero. In this case, minimizing free energy 
( |29] l with respect to A leads to constant A = 1 independent of the details of the interface, 
making A a redundant field from the point of view of any interfacial dynamics. Further- 
more, in the absence of higher derivative terms in the effective free energy expression, it 
is not possible to achieve variations in interfacial width underlining joining and splitting 
of interfaces. The interfacial tension obtained after fixing A = 1 is 7 = 2Ai. In the case 
treated in section 5, homopolymers mixture, agreement is obtained with the one proposed 
inEI. 

The proper way of deriving the equilibrium equations is by varying the metric and the 
interfacial width (namely, a{y, z)) and A(y, z)). The resulting equations are not transpar- 
ent and other types of variations, though not fully rigorous, are employed further. The 
equations below are obtained by variation of free energy in terms of A alone. The resulting 
equation describes the equilibrium conformation (minimum of the free energy) when geo- 
metric properties of the interface {H, K) are given and it is used to evaluate the interfacial 
width: 

Ai(-A-2 + 1) + 3A2A2 + 2A3X{-H) + ((A4 + A5)H^ - ^K) = (34) 

where A2..5 are correction terms as compared to Ai. The equilibrium interfacial width 
( |28T l, depends on the values of curvature invariants. For clarity, the discussion is focused 
on the interplay between the interfacial width and only one of the curvature invariants, 
namely mean curvature H. This is achieved by an additional condition, S{(j)) = in 
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equation Q, leading to ^5 = in equation ( l28T l. Variation of free energy with respect to 
both A and H gives: 

Ai (-A-2 + 1) + A^ + 2A:i\{-H) + A^H^ = (35) 
-AgA^ + 2AiXH = (36) 

Solving the equation ( [36] ) with respect with H and replacing it back in equation ( [35] l leads 
to: 

i? = -^A (37) 

-A-^ + l + 3( ^^|^^-f^ )A^ = (38) 



The exact solution is given below: 
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'4^2 ^4 - ^ 



3 



2 - 




A = 

H = 

An approximate zeroth order solution is more relevant for the following discussion. The 
coefficients accompanying the higher derivative terms are small corrections to the zeroth 
order term in the free energy. In zeroth approximation, the terms containing A2..4 lead to 
the conclusion that the constant a is a first order correction, giving: 

A = 1 

Although coefficients A3 and A^^ are first order correction terms, their ratio is not nec- 
essarily small in first order and that amounts to having a zeroth order nontrivial intrinsic 
mean curvature of interfaces. Thus, even in the zeroth order approximation, the interface 
could show a nonzero intrinsic mean curvature. In the first approximation, the inverse of 
the interfacial width departs from 1. The zeroth order term in the interfacial free energy 
density constraints A to have values of order 1 irrespective of the values of the corrective 
terms in the free energy expression In this zeroth order approximation, when A = 1, 
a connection with Helfrich interfacial free energy ([T]l can be made using equations ( [30l ) to 
(US. 

The discussion above is valid for the ideal equilibrium configuration, but a realistic dy- 
namics of field theory brings in particular points where curvature values differ significantly 
from those at equilibrium. Alternatively, the curvature is influenced by boundary condi- 
tions, insertion of high curvature boundary surfaces or external chemical potentials. For 
this reason, the first equation of (|35] l is considered such that the mean curvature H = Hq 
is fixed to a value different from that found at equilibrium and the interfacial width A is 
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obtained by solving the equation ([35]) : 

Z^X^ - 2^HoW' + (1 + ^Hl)\^ -1 = (40) 
Ai Ai Ai 

Two limiting cases are considered: flat surfaces with = Q and surfaces with high 
curvature, Hq oo: 



1 + 12^ -1 
Ho -> : A -> I .4i^ ^ I (41) 



1/2 



1/2 



Al 



1/2 



H, ^ ±^:X^(j-J (42) 

In particular, for the case when A2 <C Ai, the limit Hq — > gives A — > 1 and the 
interfacial tension of the flat interface is 2Ai. The limit H ^ 00 shows that the interfacial 
width, A~^ goes to 00, thus at high curvature points, the interface becomes smoother and 
the phases start mixing along the normal direction to the interface. Thus, high curvature- 
points are smoothed out by an increase of interfacial width to 00 (A ^ 0) leading to the 
conclusion that higher derivative terms do not allow exposed high-curvature interfaces to 
appear dynamically. 



5 A real system - The binary polymeric mixture 

A particular free energy expression that contains higher derivative terms is considered 
in this section. This expression comes from a mixture of two polymeric components, 
having the same Kuhn length a, and different polymeric lengths Na, Nb interacting with 
a coupling constant, Flory Huggins parameter, x- A. free energy of type (|2]i is derived in 
|7| for the binary incompressible mixture of homopolymers by a formal expansion of free 
energy in terms of i?g V^, where — represent either polymer radius of gyration 
(when N is replaced by particular polymeric lengths Na,b)- The formal expansion in 
powers of R^W^ is cut off at order 2, 0(2), in the expansion term and the following 
expression is obtained, see also [10|: 



kBTj dr[VM<^) + 3e^(, „ ^) + 6480 " 



«W)!/iV. JV^ (43) 
2592 V ^ 1 - ^ ' ' ^^^^ 



where fcs is the Boltzmann constant and T is the temperature. Vfh the Flory-Huggins 
potential: 

V.,„,,).*M« + (W)^l(i^ + ,44, 

The enthalpic term, x0(l ~ 0) is contained in the Flory-Huggins potential all the other 
terms being of entropic origin. Entropic terms describe two types of entropies. The trans- 
lational entropy is described by terms of Vfh depending of Na, Nb and the regular 
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(V(/))^ containing term. The entropy associated with the internal degrees of freedom and 
in particular with polymer "stretching" ability is represented by higher derivative terms in 
the free energy expression. For given polymeric lengths, a critical value for the enthalpic 
parameter XcriUcai is calculated: 




Xcriticai = 2 ( + ) (45) 



For values of x smaller than Xcriticai, the Flory-Huggins potential, Vfh has a single min- 
imum and describes a homogeneous polymeric mixture, whereas for values of x greater 
than Xcriticai, phasc segregation is possible and in fact Vfh develops two distinct min- 
ima. Another feature of the above free energy is its dependence on molecular weight of the 
polymers. As the molecular weight of homopolymers increases, there is a redistribution of 
the entropic free energy consisting in a decrease in translational entropy term, multiplied 
by — factors, and an increase in internal entropy, multiplied by Na,b factors. The en- 
tropic energy redistribution is consistent with the decrease in translational entropy as the 
molecular weight increases and an increase in internal degrees of freedom entropy as the 
molecular weight increases. 

The Flory-Huggins potential is asymmetric when A^^i ^ Nb and the two minima have 
an energy difference between them resulting in a bulk contribution to the free energy. 
The analysis of the bulk contribution is beyond the goal of this analysis and appropriate 
physical limits will be assumed to reduce the discussion to interface alone. When x is 
large enough (an order of magnitude greater than Xcriticai or roughly, x^a ^ 2 and 
xNb S> 2), the strong segregation limit is valid and although Vfh is still asymmetric, it 
presents two approximately equal minima localized at ^I^- w and ^1*+ w 1. In addition, 
since expressions from section 4 are used, the assumption of negligible variations in A is 
valid. Within such a model, the following expressions are employed to calculate yl's in 
equations (|29] l: 

V{^) = Vfh{cI>) (46) 
iV(0) = Q(0) = 5(0)=O (48) 



a2 


1 


36 0( 


1-0) 


Q(0) 


= 5(0) : 




/ Na 


6480 


[ 02 




{Na , 


2592 V ' 



The strong segregation limit leads to 5*-,+ being and correspondingly, 1 and all expres- 
sions are Taylor expanded in powers of ^^y^ ^ . The formulas developed in section 4 are 
applied to this realistic free energy and the analytic expressions are obtained using an ex- 
pansion in ^j^^ ^ . The values of and ^_ close to and 1 allow analytic integration 
over 'J for the expressions derived in this section. Calculations are done in the first two 
relevant orders for each A. Only Vfh depends on x, and the potential is rewritten in a 
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form more appropriate for perturbation: 

W*) = X*(l-«(l + ^^ + ^^^) (5.) 

V xNa (1 - cp) xNb / 

The expression is exact and it is reflected in a perturbation fashion in the following expres- 
sions: , and 4'"^, which appear in ( |29l ). These functions are treated as function of VE* 
using ([Toll, CD and gUl to (|50li: 

QJX , J 1 ln(l-*) 1 ln(*)\ 

a ^ \ 2xNb ^ ^2xNa1-^J 
1 ^ a 1/1 ln(l-^') 1 ln(*)\ 



6^ ^'(1 - *) V 2xiVs * 2xiVA 1 - ^-y 

a2 ^ X^A xNb xNa 1-* 

(l-3'J')ln(l-^')' 
^ xNb ^ 

io2 2 / 



1 1 (2 - S^-) ln(^') (1 - 3*) ln(l - *) 



XA^A xNb xNa 1-* XA^B * 

The corresponding coefficients in their first two orders in ^j^^ ^ are given below after 
evaluation of integrals ( |29] | using Mathematica l,12J and the equations above: 



1 1/2 f. TT^ 1 / 1 1 

.3 . J^„^w.,-«,,„il4„^-^„ 



A, 
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The first term OF Ai was derived [2| for a free energy having only a quadratic derivative 
term with a different numeric coefficient and without higher derivative terms. For expres- 
sions above only the first two order terms of each A are considered. For example the 
intrinsic mean curvature ( [39] l is 

i^.l^^^Xd^ (52) 
5 u^Na + Nb 

The most important parameter controlling the intrinsic curvature is the ratio between ho- 
mopolymers molecular weights: The closer this number is to 0, namely the larger 
the difference between the molecular weights the highest the value, the higher the intrin- 
sic curvature of the interface. Also, the interface curves towards the higher molecular 
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weight polymer. Similarly, there is no intrinsic curvature for a system where Na ~ Nb- 
Another conclusion is that lower values of Kuhn's parameter, a and enthalpic parameter, 
X promote higher values of intrinsic curvature. In the approximation used in this work, 
the low values achieved by x quite limited since x > XcriUcai, whereas low values 
of a are more accessible. As a decreases, the steepness of the interfacial profile 'i>{x) 
increases. Such increase, at some point, would contradict the formal perturbation expan- 
sion in i?gV'^. Lowering a will result in steeper interfaces increasing the weight of higher 
derivative terms and finally shifting the contributions from translational entropy to inter- 
nal entropy. Another conclusion regards the dependence of interfacial width of curvature. 
Dependence of interfacial width on mean curvature for large values of curvature Ho is 
dUll: 

144 1 1 
^^^xiNA + NB)\Ho\ ^^^^ 

The expression above shows that sensitivity of the interfacial width to interfacial curvature 
is increased for low values of Kuhn's parameter, a, and also for low values of x{Na + Nb) 
(within the lower limit XcriUcai)- The results highlight the intermediate regime of interest. 
Reaching the exact parameter region of interest in general form is difficult, but detailed 
parameter ranges can be obtained for particular cases. The difficulty for reaching an exact 
parameter range stems from the regime that is studied: the intermediate range where both 
translation and extension entropy play a role. The following general conclusions can be 
reached. Lower values of x, closer to XcriUcai as well as lower values of Kuhn's length 
will be more favorable for interfacial modifications. By contrast, high values for x and for 
Kuhn's length, a, decrease sensitivity increasing stability of interfaces. Larger discrepancy 
in molecular weights Na and Nb will promote interfaces with higher curvature oriented 
towards the phase with higher molecular weight. 



6 Conclusions 

In this work, high derivative free energy terms (higher than V^) for phase separated sys- 
tems (systems having domain-wall solutions) are analyzed. These terms are directly asso- 
ciated with geometric invariants, such as mean and Gaussian curvature terms in interfacial 
free energy. It is shown that these higher derivative terms allow an energetic interaction of 
geometry of the interface with the interfacial profile. In particular, it is shown that interfa- 
cial curvatures affect interfacial width, giving thus a mechanism that can be used to achieve 
in a continuous manner the joining and splitting of interfaces. The method is applied to 
a realistic free energy used in homopolymers blends and qualitative results about the pa- 
rameter range most appropriate for interfacial modifications are presented. In the case of 
homopolymer blends, higher derivative terms come from the entropy of "stretching" and 
it is shown that they affect directly the geometry of the interface. 

Finally, this work shows the potential benefits of introducing higher derivative terms in 
field free energy models used in computational modeUng of chemical systems exhibiting 
phase separation and rich membrane-like structures [T3l . The terms are potentially useful 
in tackling other types of systems where interfaces and curvature contributions appear at 
molecular scales. 
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7 Appendix 1 - Geometric free energy terms 



In this section the expressions used in section 2 are derived and the connection between 
fields living on the three dimensional space and those living on two dimensional interfacial 
space is established. Following the procedure described in section 1, the three dimensional 
field theory is written, using a field redefinition, in terms of field /. The field / has a direct 
relation with the interface (the surface where / = 0). Considering a point on the interface, 
a coordinate system can be chosen such as the interfacial point under study is (0, 0, 0) and 
a physical Ansatz for the three dimensional field / is given in section 2. The function 
is given by /(x, y, z) = A(i/, z){x — a{y, z)), where A describes the inverse of width of 
the interface at the interfacial point {y, z) = (0, 0) and a describes the location of the 
interface. Following this Ansatz, the equation for the interface in the neighborhood of 
interface is f{x, y,z) = or equivalently, x = a{y, z). Since (0, 0, 0) is an interfacial 
point, cr(0, 0) = and by an orientation choice of the coordinate system, dycr — dzcr — 0. 
The y and z coordinates are tangent to the interface at point (0, 0, 0) and the x axis is 
normal to the interface. Using the equation for the interface, the expressions of the first 
fundamental form (metric, ds^) and the second fundamental form {liap) of the interface 
at point (0,0,0) are computed using Monge representation in terms of the field a (for 
additional details see ifTTl '): 



ds' = 



(l + [dyuY^dy'^ + 2dycrdzcrdydz + (l + {d^crf^dz'^ 



(^1^° dy^ + dz' 



IIzz = = > OzzCrlo 



dyyCr 






- (5za)2 




9zzO- 








- {dzcrf 




dyzcr 






- [dyoY 4 


- {dzcjf 



TT "v^" (a,z)^o p. I 

'^'-y^ — — u ,^ =^ — ^ Oyz(y\Q 

In the same limit, (y, z) on the interface, the expressions for the curvature scalar 
invariants, mean curvature H = ^Tr{I\) and Gaussian curvature K = Det{ll), are given 
below: 

H iv^alo (54) 

K {dyya dzza ~ {dyz<7)^)\o (55) 

In order to make contact with the terms appearing in the three dimensional field theory, all 
quantities in the three dimensional field theory are computed in the limit (y, z) with 
the normal coordinate, x, being arbitrary (still in the neighborhood of the interfacial point 
(0, 0, 0)). The derivatives of the field / in this limit are shown below: 

dxf{x,y,z) = A^^-^°A|o 

dyf{x,y,z) = dy\ {x - a) - \ dya'^^'-^^ X dy\\o 
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dzf{x,y,z) 


— dzX [x - 


- cr) - A dza x dzX\o 




dxxf{x,y,z) 


= 






dyyf{x,y,z) 




— (t) — 2 dyX dya) — A dyycr ^"'—l^ ° c^j/yA 


- A Ilyy) 


dzzf{x,y,z) 


= dzzX {x 


- ct) - 2 dzX dzcr - X dzzcr (a; Q^^A - 


An,,)|o 


dxyf{x,y,z) 


^ dyX — 


^" dyX\o 




dxzf{x,y,z) 


= OzX — 


^"dzX\o 




dyzf{x,y,z) 


= dyzX {X 


— a) — dyX dzO — dzXdyCf — Xdyz<7 
{y,z)^o (a; dyzX - A ny^)|o 





where liap is the second fundamental tensor of the interfacial surface. In particular the 
Laplacian of the field / is connected with the mean curvature scalar H. 



V2/(x,y,z) 



(x V2A-2 Ai?)|c 



(56) 



Using the equations above, each of the terms in the free energy density can be connected 
to interfacial fields such as interfacial depth A and interfacial metric (namely, ct): 



(v/)2 




(A2+x2(VA)2)|o 








(X^+xH^Xf)\ 




(V/)2V2/ 




(A2+a;2(VA)2)(a;V2A- 


-2AF)|o 






xX{VXf +x'^do,XdpXd' 


'f^X- x^XdaXdpXIi°''^\Q 






{xV^X-2XH)\ 








x^do^pXd'^'^X + 2(VA)2 - 


- 2xXdc,pXIP'^ + X^li^pli'^'^ 



where the derivatives, in terms of V, are understood to act the space where each function is 
defined (i.e. if the field is /, then three dimensional V and are used, whereas if the field 
is A or CT two dimensional V and are used). In addition for the final expressions used in 
section 2, x on the right side is replaced by ^. The metric on the three dimensional space 
is flat and thus no distinction appears between upper and lower indices. The final general 
expression shows a geometric invariant character For this reason, the fully geometric 
invariant expression of each term was highlighted on the right side of the equations above. 
In addition, the term Ila^ 11"'^ can be expressed in terms of the curvature invariants H 
and K as shown in the equation (fT6] l, ( [TtI i of section 2. 



8 Appendix 2 - Dimensional reduction 

In this section, the formulas used in section 3 are derived. The free energy action in three 
dimensions is reduced in the normal direction on the interface. In order to simplify the 
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action given in the equation ( fT9t . within the general setup given in ( fTSl l. a few properties 
for the domain-wall solution of 0^ potential of section 3 are added : 

"^(x) = tanh(a;) 

= 1 - tanh(x)2 (57) 
= -2 tanh(a;) (l - t£inh{xf) 

For the first term in ( fT9] l. the term that depends on the coefficient M(^') = ^ in equation 
(fTSl l. two integrals are evaluated and the results obtained using Mathematica liT2l are listed 
below: 

/ d/*'2 = ^; / dff^' = / dff^" = ^ 

leading to equation ( |20l ). Similarly, for the second term in ( fT9] l. the term that depends on 
the coefficient N, three integrals are evaluated and the results are listed below: 

'4 _ 32, ^„2,,'4 _ 4(6^' - 49) 



35' ./ _ 315 



OO 

- 2(3^^-35^2 + 60) 

-OO 



225 



with the result given in ( |2T] ). The third term in ( fT9l ). the term that depends on the coefficient 
P, is written in (l22l l after the evaluation of seven integrals: 



df^'^-ii" = / dff^'H"f = I dff^'^^" = 

OO J —OO J — OO 

+ 00 /' + 00 

djj^'^ = / d/fvI/'3 = 

OO J — OO 



Finally, the last term from ( fT9l l. the term that depends on the coefficient R, comes with a 
contribution to the interfacial free energy given by ( l23T l after the removal of a total deriva- 
tive term and evaluation of additional six integrations: 



15' 7_oo 45 ' 225 

+00 r-\-oc p-^oo 9 

= / dff-^'-^" = 0; / d//^''*" - 



OO 
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